T-Shape Molecular Heat Pump 
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Abstract 

We report on the first molecular device of heat pump modeled by a T-shape Frenkel-Kontorova lattice. The system is a 
three-terminal device with the important feature that the heat can be pumped from the low-temperature region to the high- 
temperature region through the third terminal. The pumping action is achieved by applying a stochastic external force that 
periodically modulates the atomic temperature. The temperature, the frequency and the system size dependence of heat pump 
are briefly discussed. 
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A heat pump is a machine or device that moves heat 
from a low temperature heat source to a higher tem- 
perature heat sink by applying an external work that 
modulates the environment of the system [l|. In the 
last decade, the nanoscale heat pump attracts more and 
more attentions due to its valuable applications 0. Re- 
cent studies have carefully analyzed the physical mecha- 
nism of molecular heat pump through classical Q [3| and 
quantum || Q 0] methods. Typically, in these schemes 
a carefully-shape external force [H and a stochastic ex- 
ternal force [U periodically modulate the levels of the 
nanoobject, leading to the pumping operation. This pro- 
cess has been confirmed theoretically to be available in a 
quantum Kubo oscillator. In classical systems a Brown- 
ian noise can lead to the reverse-direction movement of 
a molecular motor Q Therefore, it is also expected 
that stochastic external force can take heat continuously 
away from the cold heat bath and pump it into the hot 
heat bath. 



a quasi one-dimensional case of this model is the forked 
nanowire and nanotube, e.g., the T-shape nanowire (l3j 
and Y-type carbon nanotube [3], which stand by much 
potential experimental execution. 
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In this letter, we build a classical prototype pumping 
device driven by stochastic external force (i.e., thermal 
noise or heat bath). In our model, as shown in Fig. la, 
the first particle of two segments, HI, 2C, connects to two 
end of segment 102 via the junction 1 and 2, respectively. 
Segment OP and 102 are coupled via the mid-particle of 
Segment 102 and the first particle of segment OP. Each 
segment is a Frenkel-Kontorova (FK) lattice. The last 
particle of segment HI and 2C are respectively connected 
to hot heat bath (high temperature, Th) and cold heat 
bath (low temperature, Tq), while the last particle of 
segment OP is the control terminal. This T-shape FK 
model, which had ever been used to study thermal tran- 
sistor and thermal logic gate [n|, is similar to the 
experimental counterpart by putting a polymer chain or a 
nanowire on the top of adsorbed sheet [12]. Furthermore, 



FIG. 1: Configuration of the heat pump. 

Segment HI (2C) and 102 are coupled to 1 (2) via a 
spring of constant k\ (fe). Segments OP and 102, which 
are coupled via the mid-particle of segment O with a 
spring of constant ka, have the same parameters with 
the exception of the size. The number of the particles 
in segment 102 is nearly twice that in segment OP. The 
stochastic external force, which is here presented by a 
heat bath, is applied to modulate the last particle of OP 
segment. The temperature of the heat bath, Tp(t), is 
a periodic square wave function. As shown in Fig.lc, 
Tp(t) equals to Tpi in the first half period and Tpi in the 
second half period, respectively. The total Hamiltonian 
of the model is 
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and the Hamiltonian of each segment can be written as 
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with XM,i and pm,i denote the displacement from equi- 
librium position and the conjugate momentum of the i th 
particle in segment M, where M stands for HI, 102, 
2C or OP. N is the number of the particles in segments 
HI, 2C and OP. The number of the particles in segment 
102 is 2N + 1. J2 H^t = Hi + H 2 + H 3 , in which Hi = 
h (x H ,i - x ,i) 2 /2, H 2 = k 2 (x ,2N+i - x c ,i) /2, and 
H3 = k 2 (xo,n+i — xp,i) 2 /2. We set the masses of all 
the particles be unit and use fixed boundaries, x\y,N+i = 
0, where W stands for H, C, or P. The main parameters 
are k m = 1.0, k 2C = 9.0, k 1Q2 = k OP = 3.0, V m = 2.0, 
V 2C = 4.5, V102 = Vop = 3.0, ki = k 2 = 0.05, k 3 = 3.0, 
T H = 0.20, and T c = 0.16 

In our simulations we use Nose-Hoover thermostat [l5| 
and integrate the equations of motion by using the 4th- 
order Runge-Kutta algorithm [16| . We have checked that 
our results do not depend on the particular thermostat 
realization (for example, Langevin thermostat). The lo- 
cal temperature is defined as Tj = (p 2 ), ( ) means time 
average. The local heat flux along the chain is defined as 
Jw,i = ku^Viixi — Xi-\)), where W stands for H, C and 
P. i is the order of the particle. The heat current, which 
flows from H to C(J H ), or O to C(J C ), or O to P{J P ), 
is defined as the positive current. The average kinetic 
energy is Ek = Yli=i~ ( v i)/(^ + 1)' where Vi is the 
velocity of the i th particle and AN + 1 is the system size, 
respectively. The simulations are performed long enough 
to allow the system to reach a steady state in which the 
local heat flux is constant along the chain. 

Here the external perturbation has the same effect as 
a heat bath coupled to the particle in the end of segment 
OP. In our model, Figure 2 shows the temperature of 
mid-particle of segment 102, Tq, changes linearly with 
the temperature of the heat bath acted on segment P. 
Due to ballistic transport of weak link systems in low 
temperature 17], the temperature of particle O nearly 
equals to the temperature of the heat bath Tp(t) i.e., 
To ~ Tp. As the temperature of the heat bath connects 
to P is low enough, the temperature of interface particle 
O (To) is also small. Therefore, when To is smaller than 
Th and To, the system absorbs energy from both heat 
baths and then the direction of the heat current is shown 
in Fig. la, here we set the heat current Jhi > 0, Jci < 0, 
and Jpi > 0. When To is larger than Th and To, the 
energy will dissipate from the system oscillator mode to 
heat baths and then the direction of the heat current is 



shown in Fig. lb, here the heat current Jr2 < 0, Jc2 > 0, 
and Jp 2 < 0. The temperature profiles along the config- 
uration of the system is shown in the inset of Fig. 2 for 
Tp = 0.02 and 0.22. Actually, the temperature of the 
controlling heat bath is variable. Provided that we use 
the appropriate values of Tp\ and Tp 2 and perform the 
simulation in long time enough (the simulation time is far 
3(2teffJ6pf)]' ian the relaxation time of the system), as shown 
in Fig.ljc, we will get a negative value of the total heat 
current in one period, i.e., Jh < and Jc < 0, which 
means a pumping operation. 
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FIG. 2: Relationship between the temperature of particle O 
and the temperature of the controlling heat bath Tp. The 
inset is the temperature profiles along segments H — 1 — O — 
2 - C as shown in Fig. la for T P = 0.02 and 0.22. 

This phenomenon can be understood from two essen- 
tial physical principles: negative differential thermal re- 
sistance (NDTR) and thermal rectification (TR) [ll| (lfl | 
17]. These two effects produce nonlinear relationship be- 
tween the heat current and the temperature difference. 
It can be explained in detail as follows: 

As shown in Figs. 3a, in the region of low temperature 
(To < To < Th), when temperature To is decreased by 
decreasing the temperature of the controlling heat bath 
Tp, Jh increases firstly and then decreases to a small 
value, however, Jc decreases linearly. The dependence of 
Jh on the temperature difference Th—To shows a NDTR 
effect. In high temperature region (To > Th > Tc), Jh 
and Jc have linear relationships with the temperature. 
In a brief, segment HO is in the open state for To > Th 
but in the close state for To < Th- This effect is not 
visible in segment OC. If we select two appropriate val- 
ues of T P , Tpi = 0.02 and T P2 = 0.22, as the temper- 
ature in the first and the second half period, the total 
heat currents Jh(Jhi + Jh2) and Jc(Jc\ + Jc2) equal 
to -4.30 x 10~ 5 and —1.15 x 10 -5 , respectively. Heat 
pump works on the condition that the negative heat flux 
is larger than the positive one during a period. 
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FIG. 3: Heat currents through two terminals H (Jh, red 
square) and C (Jc, blue circle) versus the temperature of the 
controlling heat bath Tp. (a) Numerical results for the system 
size 4N+1=65 and (b) analytical results based on Eq.(5). 



The corresponding analytical method is also included. 
As reported in Refs [lU and we replace the first 
and second derivatives of the external potential by their 
thermal average with respect to the effective harmonic 
Hamiltonian, and then equation 2 can be approximated 
by self-consistent phonon theory as 



right phonon spectra. Provided that Th and To are avail- 
able, then we can obtain the heat current flows from H 
to O, Jh- Similarly, as k 2 — > 0, we can also get the 
heat current from O to C, Jc- Due to the linear rela- 
tionship of To and Tp, figure 3b analytically confirms 
the NDTR and TR effects in T-shape FK lattices, which 
is numerically presented in Fig. 3a. 

Heat pump is a dynamical and nonequilibrium effect 
during macroscopic time. Therefore, heat pump has some 
valid conditions. We give three main parameters which 
influence the state of heat pump significantly. 

Temperature dependence In order to obtain four heat 
currents, J hi, Jh2, Jci, and Jci, shown in Figs, la and 
lb, obviously the temperature To should satisfy To < 
T c < T H in the first half period and T c < T H < T 
in the second half period, respectively. Since To almost 
equals to Tp, as illustrated in Fig. 2, Tp has to satisfy 
T P <T C < T H in the first half period and T c < T H < T P 
in the second half period, respectively. Here the main 
point is what is the range of Tp(t). Generally, we fix the 
low temperature level Tpi(= 0.02) and change the high 
temperature level Tp2. Figure 4 displays that the total 
heat currents, Jh and Jc, are negative as Tp2 ranges 
from 0.210 to 0.303. This range of Tp2 is defined as 
pumping region. When T P2 < 0.210 or T P2 > 0.303, 
pumping effect cannot be realized. It is worth mentioning 
a point, where Tp 2 = 0.251 and Jh — Jc — —0.51 x 
10 -5 < 0, indicates the optimum pumping. 
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where Tm means the average temperature and M refers 
to segments HO, OC and OP. Here we solve the tran- 
scendental equation 4 through calculating the intersec- 
tion point of left part and right part of the equation ■ 
As ki — > 0, considering classical Landauer-type equa- 
tion, we can get the heat current flows from H to O as 
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in which the transmission coefficient is 



(5) 



FIG. 4: Temperature dependence of the heat current in a T- 
shape FK lattice. The vibration period of Tp(t) is 5.0 x 10 6 
steps and the system size is 65. 
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The cutoff frequencies range from ui\ 

to U>2 
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Frequency dependence In our simulation, the system 
relaxation time (r s ) of our model is about 6.0 x 10 9 simu- 
lation steps, then r s =6.0 x 10 . The oscillating frequency 
of single particle (ui p , or l/r p ) ranges from 2.5 x 10 -3 to 
6.2 x 10 _1 , which can be calculated through the phonon 
spectral analysis of the particle's velocity [10(. In order to 
get a stable heat current, the vibration period of temper- 
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ature level Tp(t) should be far smaller than the system re- 
laxation time. Furthermore, the vibration period of noise 
level cannot be near the oscillating period of single parti- 
cle. As shown in Fig. 5a, in the case of low frequency, the 
heat pump works normally with a stable negative heat 
current Jh- However, when the vibration frequency of 
temperature is increased to a value larger than 3.0 x 1CT 3 , 
the direction of heat current changes and the heat pump 
stops working. It is easily understood from the change of 
the average kinetic energy with the frequency. As shown 
in Fig. 5b, the system maintains its average kinetic en- 
ergy onto a value E K (= \{E{T P1 ) + E(T P2 )) = 0.072) 
at low frequency. However, the average kinetic energy 
has a significant change in the frequency region of sin- 
gle particle, which corresponds to the shadow area in 
Fig. 5. In this region, the frequency may match the os- 
cillating frequency of single particle, the temperature of 
the controlling heat bath has a significant influence on 
the oscillating energy of every particle and the system is 
nonconservative. Therefore, we set the vibration period 
of temperature level, r„, be 5.0 x 10 4 , which satisfies 

Tp << T„ << T„. 



and then NDTR disappears. Thus, the valid condition 
for heat pump will be unavailable. 
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FIG. 6: The finite size effect of the pumping efficiency for 
two kinds of heat baths, Langevin (solid Squares and Circles) 
and Nose-Hoover (open Up Triangles and Down Triangles). 
The remaining parameters are Tp\ = 0.02, Tp2 = 0.22, and 
r n = 5.0 x 10 6 . 
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FIG. 5: Frequency dependence of (a) heat currents and (b) 
average kinetic energy of heat pump. The remaining param- 
eters are the same as for Fig.4. 

System size effect The results displayed above are for 
a system with 65(4N+1) particles. Since the heat pump 
mechanism in our model is due to the coupling between 
two asymmetric lattices it is reasonable to expect that 
the system size will definitely influence the heat pump 
efficiency. As shown in Fig. 6, the heat current is de- 
pendent on the system size. The heat current of heat 
pump increases firstly and then decreases by increasing 
the number of particles. Finally, the heat pump stops 
working when the system size is larger than 500. This 
phenomenon can be understood by system size depen- 
dence of NDTR [ItJ ■ When the system size is increased, 
the system goes to completely diffusive transport regime 



We have to point out that in this letter we use the more 
complex four segments rather than three segments FK 
lattices just for the convenience of theoretical analysis. 
Actually, in the case of three segment lattices, the heat 
pump may work more efficiently. Finally, we would like 
to discuss the improvement of our heat pump model. Fig- 
ure 1 displays a single heat pump working between two 
heat baths with small temperature difference. Moreover, 
this kind of single heat pump works with low pumping 
efficiency. Therefore, if we expect to get a more powerful 
heat pump, we can connect single heat pumps in series 
or in parallel. 

Up to now we only pay attention to the behavior of 
the device as a heat pump, which requires Jh < and 
Jc < 0; however, for the situation with only Jc < 
(when the system size is larger than 500 as shown in 
Fig. 6), which would be useful as a refrigerator mode to 
extract heat from the cold source, although in this case 
this heat will not go to the hot source but will go to the 
oscillating heat bath; anyway, the ensemble will globally 
act as a refrigerator. 

In conclusions, we have reported the feasibility to pro- 
duce molecular heat pump based on T-shape FK model. 
This device can pump the heat from the low-temperature 
region to the high-temperature region through control- 
ling the atomic temperature of the third terminal. Al- 
though the heat pump presented here is only an ideal 
model, it can be easily imitated in experiment. The 
study may also be a valuable illumination in fabricating 
nanoscale heat pump. Besides, our heat pump model will 
help deeply understand the effect of negative differential 
thermal resistance. 
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